CLAIMS 

1.-52. (CANCEL) 

53. (NEW) A method of estimating arterial delay and arterial dispersion (t, a, a) 
values for outputting blood perfusion indices for a region of interest (ROI) by 
operating a computer program on intensity data in a computer comprising: 

a. applying a first gamma-variate function (GVF) to an arterial input function 
(AIFa) to provide an estimated first model of a vascular transport function 

ha(t), wherein for t < ti, ha(t) = 0 and for t > t,, ha(t) =—{t- 1^' e~^'~''^^''' , 

Ox 

wherein an estimated ti is the transit time of a contrast agent from a measured 
initial said AIFa to a region of interest (ROI); 

b. estimating an initial value Oi of said contrast agent, wherein said 
ai=(ti)(Pi)/(l-Pi), wherein said Pi is a known relative dispersion value having 
a range from 0 to 1; 

c. convolving AIFa(t) with said ha(t, ai=0) for obtaining an arterial input 
function AIFt(t) = AIFa(t) ® ha(t, ai=0) at said ROI; 

d. estimating a blood flow rate Ft and a tissue impulse residue function R€(t) by 
deconvolving a concentration curve C(t) = (Ft/kH)AIFt(t) ® KJfy, wherein kn 
is a hermocrit correction constant having a known value; and 

e. outputting estimated and optimized tissue mean transit time and dispersion (t2, 
CC2, Oi) values from an estimated transport function he(t) for input to a 
simulated transport function hs(t), wherein a simulated tissue impulse residue 
function Rs(t) is determined, wherein a simulated concentration curve Cs(t) is 
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fitted to said measured C(t) and quantitative said blood perfusion indices are 
calculated. 

54. (NEW) The method of claim 53, wherein said intensity data is generated by 
administering a contrast agent to a body lumen of a body during a dynamic 
imaging scan, wherein said body lumen comprises an artery or vein, wherein 
an image response from said contrast agent is recorded to computer data 
storage in a computer. 

55. (NEW) The method of claim 53, wherein said C(t) is a temporal concentration 
of said contrast agent obtained from said intensity data, wherein said intensity 
data comprises contrast images sequentially acquired from a region in a body, 
whereby said contrast agent concentration is plotted versus time. 

56. (NEW) The method of claim 53, wherein said AIFa is based on a measured 
early arrival contrast agent peak intensity from a feeding blood vessel to said 
ROL 

57. (NEW) The method of claim 53, wherein said AIFa is scaled upward according 
to a venous input ftmction (VIF), wherein said VIF is based on a measured late 
arrival contrast agent peak intensity from a large vein draining from said ROL 

58. (NEW) The method of claim 53, wherein said estimated transit time ti is the 
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transit time of said contrast agent from a measured initial said AIFa of said 
contrast agent C(t) in a body lumen to said ROI, wherein said ti is estimated 
from plots of said AIFa versus time and said C(t) versus time. 



5 59. (NEW) The method of claim 53, wherein said ha(t) is calculated using said 

estimated transit time ti and said estimated dispersion value Oi, wherein ha(t, 
ai=0) is plotted versus time. 



60. (NEW) The method of claim 53, wherein said estimated transport fimction he(t) 
10 is calculated using the relation he(t) = - dRe(t)/dt. 



61. (NEW) The method of claim 53, wherein said tissue mean transit time and 
dispersion (t2, a2, 02) values are estimated from said estimated transport 
fimction he(t), wherein said ta, said 02 and said a2 are input to a simulated 

15 transport fimction hs(t), wherein said hs(t) is said second gamma-variate 

function. 

62. (NEW) The method of claim 53, wherein said simulated tissue impulse 

resistive fimction Rs(t) is determined using the relation Rs(t) = 1 - ^hs{x)dx, 

0 

20 

63. (NEW) The method of claim 53, wherein said simulated concentration curve 
Cs(t) is determined using the relation Cs(t) = (Ft/kH)AIFt(t) ® Re(t) = (Ft/kn) 
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t 

/ AIFt(t) Rt(t-T)dt. 

0 



64. (NEW) The method of claim 53, wherein said Ft, said ti, said d, said a2, said 
t2, said a2, and said 02 are optimized by a least squares method to fit said Cs(t) 
to said C(t). 

65. (NEW) The method of claim 53, wherein said perfusion indices have the 
relations: 

a. blood flow (BF) = Ft; 

b. Mean Transit Time (MTT) = t2 + 02(1 ■i-a2); 

c. Blood Volume (BV) = BF* MTT; 

d. Arterial Delay (DT) = ti + ai(l+ai); 

e. Arterial Dispersion time (ADT) = 0\ Vl + oci ; 

f. Tissue Dispersion Time (TDT) = 02 Vl + aa ; 

g. Relative Arterial Dispersion (RAD) = ADT/DT; and 

h. Relative Tissue Dispersion (RTD) = TDT/MTT. 

66. (NEW) The method of claim 53, wherein said AIFt(t) is measureable in a small 
lumen showing a delay relative to said AIFa(t), wherein optimized values for 
said oi and said ti are determined by fitting said simulated AIFt(t) to said 
measured AIFt(t), wherein said relative dispersion Pi is determined and appUed 
to all other said intensity data of said ROI using said Pi, wherein a more robust 
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fitting process is provided by a reduced number of parameters for optimization. 

67. (NEW) The method of claim 66, wherein when said relative dispersion Pi 
is determined, said vascular transport function ha(t) is described by a single 
variable said ti with a constant said pi, wherein a two-step method is used 
to determine said delay and said dispersion values comprising: 

a. deriving an initial tissue impulse residue function Ro(t) by 
deconvolving C(t) = (Fo/kH)AIFa(t) ® Ro(t) using a model-firee singular 
value decomposition (SVD) method, wherein said time delay ti is 
determined by a maximum position of said Ro(t) at Ro max(t==ti); and 

b. determine said AIFt(t) at an input of said ROI using said ha(t) with said 
ti and said Pi held constant, wherein said 0\ is determined. 

68. (NEW) The method of claim 67, wherein a value of tissue blood flow Ft 
15 and a corrected impulse residue function Re(t) are obtained by 

deconvolving C(t) = (Ft/kH)AIFt(t) ® Re(t) using said SVD method, 
wherein said perfusion indices are determined from a curve of said 

Re(t), wherein MTT= /J Re(T)dx, BF=Ft, and BV=BF*MTT. 

20 69. (NEW) The method of claim 53, wherein said contrast agent is in a tissue ROI 

having a tissue mean transit time r, wherein a tissue impulse residue function is 
approximated by the relation R(t >r) = Ee"*^^^"^^ and R(t^T) = 1, wherein E is an 
extraction fi-action of said contrast agent in said tissue, wherein k is a constant 



5 
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clearance rate of said contrast agent diffusing from said tissue having a relation 
k = E*FtA^e, wherein Vcis the volume fraction of extravascular and 
extracellular space (EES) in said tissue. 

70. (NEW) The method of claim 69, wherein said tissue impulse residue 
fimction Rs(t) of said simulated concentration curve Cs(t) is replaced by an 
average impulse residue function that incorporates said contrast agent 
leaked out of a blood vessel into said tissue and gradually clearing from 
said tissue, wherein said simulated concentration curve Cs(t) is fitted to said 
measured C(t) and quantitative said blood perfusion indices are calculated, 
wherein said E and said Ve are additional parameters optimized with other 
adjustable parameters using a least squares method. 
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